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Using the average method, we derive a closed set of linear 
equations that describes the spinning up of an harmonically 
trapped gas by a rotating anisotropy. We find explicit expres- 
sions for the time needed to transfer angular momentum as 
well as the decay time induced by a static residual anisotropy. 
These different time scales are compared with the measured 
nucleation time and lifetime of vortices by the ENS group Q . 
We find a good agreement that may emphasize the role played 
by the non-condensed component in those experiments. 



condensed component. Finally, in Sec. VI, we consider a 
related problem: what is the time needed to dissipate a 
given angular momentum by a static residual anisotropy 
? We have in mind the role of the axial asymmetry of a 
magnetic trap (Ioffe-Pritchard or time-orbiting potential 
traps) induced by the presence of gravity. We show that 
this effect may explain the finite lifetime of vortices. 



I. A REMINDER ON THE ROTATING FRAME 



INTRODUCTION 

Superfluidity of a Bosc-Einstcin condensate is natu- 
rally investigated by its rotational properties 0J§]. So 
far, two experimental schemes have been successfully ap- 
plied to gaseous condensates in order to generate quan- 
tized vortices . The first one uses phase engineering 
by means of laser beams whereas the second one is the 
analog of the "rotatingbucket" experiment, initially sug- 
gested by Stringari [||J6|. 

In the latter, atoms are first confined in a static, axi- 
ally symmetric Ioffe-Pritchard magnetic trap upon which 
a non-axisymmetric attractive dipole potential is super- 
imposed by means of a strirring laser beam. In this pa- 
per, we address the question of the behaviour of the non- 
condensed component for this experiment. We investi- 
gate the transfer of angular momentum to an ultracold 
harmonically confined gas by such a time dependent po- 
tential. We also study possible mechanisms for its dissi- 
pation. First (Sec. I), we recall the Lagragian and the 
Hamiltonian formalism for a single particle in a rotating 
frame. Then in Sec. II, we give the expression for the 
rotating potential that has been used in our study. We 
briefly expose analytical results for the single particle tra- 
jectory in this potential in Sec. III. The rest of the paper 
deals with the crucial role played by collisions. A classi- 
cal gas that evolves in this potential thermalizes in the 
rotating frame, leading to a finite value of its mean an- 
gular momentum. We investigate this equilibrium state 
in Sec. IV. In Sec. V, we derive the analytical expression 
for the time needed to spin up a classical gas with an ap- 
proach based on the average method [Q . We evaluate the 
time needed to transfer angular momentum and thus we 
deduce the characteristic time for vortex nucleation via 
angular momentum transfer from the uncondensed to the 



In this section, we recall the hamiltonian for a classical 
particle in a rotating frame characterized by the fixed 
rotation vector ft. Without loss of generality, we choose 
a rotation around the z axis ft = fle z and the same origin 
for the rotating frame Sft' as the one of the laboratory 
frame 5R (see Fig. 1). In the following, quantities with a 
prime are evaluated in the rotating frame. Coordinates 
in 5ft' are linked to coordinates in 5ft just by a rotation of 
angle 8 — fit: 
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The lagrangian for the single particle movement in the 
rotating frame reads M: 
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m(fl x r) 2 + mv'.(f2 x r) - V el 
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where V ex ± is the potential energy that describes the role 
of external forces. The correspondance between the lab- 
oratory and rotating frames for momentum, angular mo- 
mentum and the hamiltonian is given by: 



, dL 1 

L' = rxp' 
W = p'.v' - 



mv + mfl x r 

=rxp=L 
£' = H-fl.L 



(3) 



Note that the angular momentum as well as the momen- 
tum are the same in both frames, but the link between 
the momentum and the velocity differs. Using the hamil- 
tonian formalism, one can easily extract the equation of 
motion in the rotating frame: 



in- 



dt 



(4) 



where F cxt refers to the external force derived from the 
potential V c ^t, F ccn is the centrifugal force: 
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F ccn = -mSl x (n x r) = - W cen (5) 

with V^ cn = — mVL 2 (x 2 + y 2 )/2, and F cor is the Coriolis 
force: 

F cor = -2mfl x v'. (6) 

Note that all formulas of the system (|J) are still valid 
even for a time-dependent rotation vector O. 



We find a window of instability around the value Q = loq, 
i.e. for £1 6 [woVl — e; uj \/1 + e]. For < woVl — e i 
the stability is essentially ensured by the harmonic trap- 
ping even if reduced by the centrifugal force, whereas for 
il > + e the Coriolis force plays a crucial role in 

the stabilization of the trajectory. The latter effect is 
similar to magnetron stabilization in a Penning trap for 
ions ffj. 



II. THE ROTATING TRAP 

In Bose-Einstein condensation experiments, the mag- 
netic confinement is axially symmetric. In order to spin 
up the system, one possibility consists of breaking this 
symmetry by superimposing a small rotating anisotropy 
as initially suggested by S. Stringari ||. This breaking of 
the rotational invariance of the external potential can be 
carried out experimentally by adding a rotating stirring 
dipolar beam to the magnetic field of the trap Q. This 
combination of light and magnetic trapping induces the 
following harmonic external potential (see Fig. 1) M: 
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{\ + e)x 12 + {l-e)y 12 + \ 2 z 2 



(7) 



where we have defined the geometric parameter A = 
Wz/wo, which is responsible for the shape of the cloud. 
For instance, if we set e = 0, the trap is isotropic for 
A = 1, cigar-shaped for A <C 1, and disk-shaped for 
A ^> 1. The potential ([?]) is time-dependent if expressed 
with laboratory coordinates (x, y, z), but static in the 
rotating frame, i.e. in terms of (x' , y' , z). 



III. SINGLE PARTICLE TRAJECTORY 

Let us first investigate the single particle trajectory. 
Equations in the laboratory frame and for the transverse 
coordinates are given by: 



cos2fi< sin2£lt 
sin2f2i — cos2fW 



(8) 



In the rotating frame, the corresponding equations are 
time independent. They can be rewritten by means of 
the complex quantity: £' = x' + iy': 



£' + 2M£' + {uj 2 - n 2 )? + ewgf = 



(9) 



where £ denotes the complex conjugate of £. The sec- 
ond term in Eq. (||) accounts for the Coriolis force, the 
centrifugal force leads to a reduction of the harmonic 
strength (third term), and the last term is the contribu- 
tion of the asymmetry. The stability of the single particle 
movement is extracted from the dispersion relation of (JTJ). 



IV. EQUILIBRIUM STATE OF THE GAS 

In practice, experiments are carried out in the so-called 
collisionless regime, since on average an atom undergoes 
less than a collision during a transverse oscillation period. 
However, collisions are of course essential to explain the 
dynamics of the gas induced by the potential (J?j). We 
consider the situation in which a gas at a given temper- 
ature To is initially at rest in the lab frame, and at t = 
the axial asymmetry e is spinned up at a constant angu- 
lar velocity fi as explained above. If < — e one 
expects that elastic collisions will ensure the thermaliza- 
tion of the gas in the rotating frame. This equilibrium 
state is defined since a minimum of the effective potential 
V e s — V cxt + V ccn always exists in this range of values for 
il. In Fig. 2, we compare the stability of a single particle 
(upper part of the diagram) with that of the interacting 
gas. In the rotating frame, one can compute equilibrium 
quantities by means of the Gibbs distribution p which 
reads Ol: 



p(r / ,v / )oce- w '< r '' v '>/fc3T q 
where TC is given by: 
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As regards statistical properties of the gas, the rotation 
is equivalent to a reduction of the effective transverse 
frequencies of the trap due to the contribution of the 
centrifugal force. The Coriolis force plays no role for the 
equilibrium state. During the thermalization, the mean 
angular momentum per particle increases from zero to 
its asymptotic value (L z ). This last quantity is easily 
derived from the Gibbs distribution (|l0|) : 

Ni ~ n ) - e w o 
For typical experimental parameters loq/2t: = 200 Hz, 
Q = ujq/2, e = 0.05 and To = 1 /xK, the angular momen- 
tum per particle (L z )/H ~ 150 unlike the superfluid part 
for which for instance the angular momentum is equal to 
h per particle in the presence of one vortex. 

To check that the gas undergoes a full rotation, one 
can search for a displacement of the critical temperature 
induced by centrifugal forces ||, or perform a time-of- 
flight measurement. In the latter case, one expects that 
the ratio between x and z size scales as (1 — SI 2 /loq)^ 1 ^ 2 
for long time expansion. 
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V. TIME NEEDED TO REACH EQUILIBRIUM 

In this section we derive the expression for the time t up 
needed to reach the equilibrium state in the rotating co- 
ordinate system. In other words, t up corresponds to the 
time needed to build correlations between x, v y , y and 
v x . To estimate this time, we use the classical Boltzmann 
equation. Our analysis relies on the use of the average 
method as explained in H. For instance, the equation 
for (x'Vyi — y'v x i) involves (x'y'} which itself is coupled 
to (x'vyi + y'v x i) and so on. Terms that do not corre- 
spond to a conserved quantity in a binary clastic collision 
lead to a non zero contribution of the collisional integral. 
For instance, a mean value such as (v x >v y i) involves the 
occurence of quadrupole deformations in the velocity dis- 
tribution which make the contribution (v x /v y i I co ii) 7^ 0. 
Here, I co \i stands for the collisional kernel of the Boltz- 
mann equation At this stage, we perform a gaussian 
ansatz for the distribution function. After linearization 
this quadrupolar contribution results in the so-called re- 
laxation time approximation H|l~3J] : 



(v x >v y >I C ou) = - 



\V X 'Vyl 



time is very long t^p ~ 15 s. Note that to transfer just 
h of angular momentum per particle, one needs only 100 
ms. This time is on the order of the nucleation time 
for vortices that has been experimentally observed by 
the ENS group (w| . One should nevertheless be carefull 
since the non-condensed component is actually a Bose gas 
rather than a classical one that evolves in a non harmonic 
potential because of the mean field potential due to the 
condensed component. 

In the hydrodynamic limit, the characteristic time for 
spinning up is: 

C = l/(2 e 2 „ 2 r). 

So far this regime is not accessible for ultracold atom ex- 
periments since inelastic collisions prevent the formation 
of very high density samples. We recover here the spe- 
cial feature of the hydrodynamic regime, i.e. the time 
needed to reach equilibrium increases with the collisional 
rate. In Fig. 4 we have reported the evolution of t up 
as a function of loqt from a numerical integration of the 
13 x 13 system. The smallest value is obtained between 
the collisionless and hydrodynamic regimes, as is usual 
for the relaxation of a thermal gas M . 



This method leads to a closed set of 13 x 13 linear equa- 
tions (see Appendix), and provides furthermore an ex- 
plicit link between the relaxation time r and the colli- 
sional rate in the sample. 

It is worth emphasizing that the dynamic transfer of 
angular momentum to a classical gas by rotating a su- 
perimposed axial anisotropy involves a coupling between 
all quadrupole modes. The average method is fruitful in 
the sense that it yields a closed set of equations when one 
deals with only quadratic moments, namely: monopole 
mode, scissors mode, quadrupole modes, ... Non-inertial 
forces are linear in position or velocity, and thus give rise 
only to quadratic moments using the average method. 

Although unimportant for equilibrium properties, the 
Coriolis force plays a crucial role for reaching equilib- 
rium. Fig. 3 depicts a typical thermalization of the gas 
in the rotating frame, leading to the equilibrium value 
( |l2| ) of the angular momentum, obtained by a numeri- 
cal integration of the 13 x 13 set of equations. After a 
tedious but straightforward expansion of the dispersion 
relation, one extracts the smallest eigenvalue that drives 
the relaxation for a given collisional regime. First we fo- 
cus on the regime in which experiments are performed: 
collisionless (ujqt ~ 10) and with e 2 <§C f2 2 /w 2 . Denoting 
t^p as the time needed to spin up the gas in this regime, 
we find: 



.CL 



(13) 



The result is independent of the geometrical aspect ratio 
of the trap A as physically expected. Using the same 
numerical values as in the previous, we deduce that this 



VI. TIME NEEDED TO DISSIPATE A GIVEN 
ANGULAR MOMENTUM 

Consider a gas with a given angular momentum, ob- 
tained for instance as explained before. If this gas evolves 
in an axially symmetric trap, the angular momentum is 
a conserved quantity. In contrast, if a small asymmetry e 
exists between the x and y spring constants, the angular 
momentum is no longer a conserved quantity and it is 
thus dissipated. We call tdown the typical time for the 
relaxation of the angular momentum. This problem is 
very different from the one we faced previously since the 
rotating frame is not an inertial frame. We thus expect 
tdown 7^ t up . As in the previous treatment, this problem 
can be computed with the average method. The cor- 
responding equations are nothing but the scissors mode 
equations ||, i.e. a linear set of 4 equations involving 
(xy), (xvy-yv x ), {xv y +yv x ) and (v x v y ) (see Appendix) . 
Searching a solution of this system of the form exp(— At), 
one finds 
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47 



1- v/l-e 2 A 2 



(14) 



where the critical anisotropy is related to r and loq by 
e c = 1/(4:UJ t). In Fig. (5), we plot different curves for 
the relaxation of angular momentum depending on the 
value of e with respect to e c . For e < e c (long dashed 
line), one has a purely damped relaxation; in the limit- 
ing case e <C e c , idown — l/(2e 2 Wgr). On the contrary 
for e > e c (solid line), one has a damped oscillating be- 
haviour; in the limiting case e 3> e c , idown — 4r, and the 
oscillating frequency is equal to eujQ. Moreover, as we use 
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a linear analysis, this decay time does not depend on the 
specific initial value of the angular momentum. For the 
experiment described in M , the axial asymmetry induced 
by gravity is on the order of 1% and the collisional rate 
is such that ujqt ~ 10. The corresponding decay time of 
the angular momentum is then evaluated from Eq. (|l4) 
to be ~ 500 ms. This value matches surprisingly well the 
typical observed lifetime of a vortex. Therefore, one pos- 
sible interpretation may be that the thermal part acts as 
a reservoir of angular momentum and thus ensures the 
stabilization of the vortex as long as this thermal part is 
itself rotating significantly. On the opposite when the an- 
gular momentum of the thermal part is zero, the vortices 
disappear. 



CONCLUSION 

In this paper we have investigated the dynamic of 
transfer and dissipation of angular momentum for a clas- 
sical gas and in all collisional regimes. We derive very 
different time scales for both processes. The short time 
deduced to dissipate angular momentum suggests a high 
sensitivity of the experiment described in W to residual 
static anisotropy (see different time scales between Figs. 
(3) and (5)). In practice, as underlined above one can- 
not avoid residual fixed anisotropy due for instance to 
gravity. The competition between a rotating and a static 
anisotropy may explain why no evidence for a full rota- 
tion of the classical gas was experimentally reported in 
Ref 1. 



ACKNOLEDGEMENTS 

I acknowledge fruitful discussions with V. Bretin, F. 
Chevy, J. Dalibard, K. Madison, A. Recati, S. Stringari 
and F. Zambelli. I am grateful to Trento's BEC group 
where part of this work was carried out. 




FIG. 1. Rotating frame with respect to laboratory frame. 
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FIG. 2. Stability diagram for the single particle and for an 
interacting gas in the potential with a rotating anisotropy. 
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FIG. 3. Mean angular momentum per particle as a function 
of time. 
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FIG. 4. Time needed to reach equilibrium by rotation of a 
small anisotropy. A minimum of t up is obtained in-between 
the collisionless and the hydrodynamic regimes. 
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FIG. 5. Decay of angular momentum due to a small resid- 
ual anisotropy. We distinguish three different regimes: e > e c 
damped oscillations (solid line), e — e c (small dashed line) the 
frontier with the region of purely damped decay (e < e c , long 
dashed line). 



APPENDIX A: APPENDIX 

Hereafter, the 13x13 closed set of equations that de- 
scribes the dynamic of a classical gas induced by the ro- 
tating potential (0): 



d(x' 2 + y' 2 + z' 2 ) 
dt 

d{x' 2 + y' 2 - 2z' 2 } 
dt 



d(x'y') 
dt 
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d{x'v y < - y'v x i) 
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2u 2 e(x'y' 
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+ — {x v x > + y iy + z v z 
+^-(x'v x > + y'v y f - 2z'v z >) = 



d(x'v y i + y'v x >) 
dt 



+ 2(u% - n 2 )(x'y' 



+2Q(x'v x > - y'v y t) - 2(v x iv y >) = i 
d{x'v x , - y'vyl + _ q2){x , 2 _ yl2 



dt 
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+ + y» + Z V) _ 2 Q{x'v yl + y'v x 

+ ^(x' 2 +y' 2 -2z' 2 )-(v 2 x ,-v 2 y ,) = 



d(x'v x > + y'vyi + z'v x t) 



, 2/12 12 

+ euj {x -y 
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2 - 2tt 2 + X 2 



i - n 2 - a 2 



lu 2 (x' 2 +y' 2 + z'- 
ul{x 12 +y /2 -2z r 



-{v 2 x , + v 2 y , + vl) - 2fl(x\, - y'v x ,) 
d(x'v x > +y'v yl - 2z'v z <) 



i 2 1 12 12 

+ ecv {x -y 
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+ f^0(l_A 2 -n 2 )(x' 2 +y' 2 + z' 2 
3 

+ ^(l + 2A 2 -r! 2 )(x' 2 +t/' 2 -2z' 2 
-{v 2 x , + v 2 y , - 2v 2 z ,) - 2n(x'v y , - y'v x ,) = I 

+ euKx'vy, - y>v x ,) + 2n(v 2 x , - v 2 y , 



(A5) 



(A6) 
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(A9) 



+ (co 2 - n 2 )(x'v y/ + y'v x ,) = _i^V) (A10) 
d(v 2 x ,-v 2 y ) 



dt 



+ 2(uj 2 - fl 2 )(xv x < - y'v y ') 



-8Q(v x 'V y >) + -^-{x'v x > + y'vyi + z'v z >) + 



^-(x'v x , + y'vy, - 2z'v z ,) = - {Vx ' Vy>) (All) 
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d(v 2 x , + v 2 y , + yp 
dt 



+ 2eu>l(x'v x > - y'vyi) 



2uj 2 



+ _£ (2 + A 2 - 2n 2 )(x'v x , + y'vy, + z'v z ,) + 
2uj 2 

-^(1 - A 2 - n 2 )(x'v x , + y'vy, - 2z'v z ,) = (A12) 



d(v 2 x , + v 2 y - 2v 2 ,} 
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+ 2eu)l(x'v x > - y'v^) 
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(A13) 



One can show from the gaussian ansatz that the relax- 
ation time t is the same for all quadrupolar contributions. 
In this system, linear terms in account for the Coriolis 
force whereas quadratic terms in f2 are due to centrifugal 
force contributions. All quadrupolar modes are involved 
in this system. In order to enlight the physics of this 
system, let us consider some limiting cases. One can 
check the conservation of energy in the rotating frame 
d(H')/dt = 0. The stationary state, obtained by setting 
all time derivatives of moments to zero, is nothing but 
the equipartition law (see Q, §44): (1 — S7 2 + e){x' 2 ) = 
(l-Q 2 - e)(y' 2 ) = X 2 (z 2 ) = (v 2 x ,) r {v 2 yl ) = (v 2 ,). If 
e = and Q = 0, the trap is axially symmetric and 
one recovers the conservation of the angular momen- 
tum. In addition the system gives rise to 3 independent 
linear systems on quadrupolar quantities: one for the 



= 2, (xy) mode, another one for m = 2, (x 2 



y 



and finally the m = mode that describes the cou- 
pling between monopole ((r 2 )) and quadrupole mode ( 
to = 0, (x 2 + y 2 — 2z 2 )) M. Finally, one can recover the 
scissors mode equations for a classical gas by considering 
the case Q = and e f p|: eqs Al, A5, A6, A10. 
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